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ABSTRACT 

Accelerating superluminal motions were detected recently by multi-epoch 
Very Long Baseline Interferometry (VLBI) observations. Here, a Linear Acceler- 
ating Superluminal Motion (LASM) model is proposed to interpret the observed 
phenomena. The model provides a direct and accurate way to estimate the view- 
ing angle of a relativistic jet. It also predicts that both Doppler boosting and 
deboosting effects may take place in an accelerating forward jet. The LASM 
model is applied to the data of the quasar 3C 273, and the initial velocity, accel- 
eration and viewing angle of its three components are derived through model fits. 
The variations of the viewing angle suggest that a supermassive black hole binary 
system may exist in the center of 3C273. The gap between the inner and outer 
jet in some radio loud AGNs my be explained in terms of Doppler deboosting 
effects when the components accelerate to ultra-relativistic speed. 

Subject headings: galaxies: jets — techniques: interferometric — quasars: indi- 
vidual (3C 273) 

1. Introduction 

Superluminal motions were first detected by Very Long Baseline Interferometry (VLBI) 
in 1970s (Cohen et al. 1977). The observations carried out since then showed that such mo- 
tions are quite common in blazars. Many models were proposed to explain the phenomenon, 
which were reviewed by Blandford et al. (1977). The relativistic jet model (Blandford & 
Konigl 1979) has become the de facto standard model to describe the superluminal motion. 
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Prom a simple physical consideration, acceleration is required to drive the subrelativistic 
material in the accretion disk up to the relativistic speed observed in the jets. Acceleration 
is also necessary to compensate the radiation losses in jets. There is increasing evidence that 
the central engine of acceleration in AGNs is a rotating supermassive black hole surrounded 
by a geometrically thin accretion disk, which gives rise to the formation of a pair of relativistic 
jets. The mechanism of acceleration may be due to the centrifugal and shear effects (Webb 
et al. 1994; Rieger & Mannheim 2002). In recent years, multi-epoch VLBI monitoring 
indeed led to the discovery parsec-scale accelerating superluminal motion in the quasar 3C 
273 (Krichbaum et al. 2001) and in other sources (Homan et al. 2001). 

Currently, quadratic or cubic polynomials are usually employed to fit the accelerating 
superluminal motion curves (Krichbaum et al. 2001). Such mathematical fits, however, have 
no real physical meaning. Actually, a physical model is required to describe the observed 
phenomena. Here, we consider a simple model where a component is ejected out coUinearly, 
with intrinsic constant acceleration g. Such acceleration may exist in rotating magnetized 
jets (Rieger & Mannheim 2002), especially at the outsets of jets where the resistance of the 
intergalactic medium (IGM) can be neglected. Also, some theoretical work shows that the 
jets in pulsars indeed move outward with linear acceleration (Contopoulos & Kazanas 2002). 
In the frame of special relativity, this model can interpret the accelerating superluminal 
motions, and we name it as a Linear Accelerating Superluminal Motion (LASM) model. 

In this Letter, the LASM model will be introduced in detail. The basic kinematic 
equations of an accelerating relativistic component and their solutions arc listed in Section 
2. In Section 3, the proper motion of 3C273 monitored with VLBI is model-fitted, and the 
viewing angle, initial velocity and acceleration of three components are estimated. Possible 
application of the LASM model, with emphasis on the estimates of the viewing angle and 
the evolution of the Doppler factor, are discussed in Section 4. 

Throughout this Letter, we adopt a fiat, accelerating cosmology with Hubble constant 
Ho = 71km s-^ Mpc-\ Oa = 0.73 and = 0.27. 

2. Linear Accelerating Kinematics 

2.1. Basic Equations 

We follow the relativistic jet model (Blandford & Konigl 1979). In this model, a com- 
ponent moves out along a coUimated jet forming a viewing angle 9 with the line of sight. 
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The kinematics of this component can be described by the following equations 

d/J-oh — sin^ dr (1) 



dt 
dtob 
"df 



9 (3) 
l-/?cos^ (4) 



/5app = (5) 
"Job 

where r is the displacement between the component and its origin, fi^h is the observed 
proper motion, t and tob ai'e the intrinsic time and the observer's time respectively , f3 and 



7 = 1/y 1 — denote the component's velocity (in units of the speed of light) and the 
Lorenz factor, g is the acceleration, and /?app is the apparent velocity. 

In order to interpret the accelerating superluminal motions observed by VLBl, the 
relationship between ii^h and tob should be derived based on the above equations. In case of 
hnear acceleration where 9 remains constant, the solution can be deduced through the solving 
routines of normal differential equations. It is worth noting that parametric equations, which 
describe the relationship between //ob and toh, are simpler to be derived than the "direct" 
solution of /iob in terms of toh- 



2.2. Solutions 



Assuming that g is constant, parametric solutions of Hoh and tob can be derived. They 



are 



sin 6 



fJ'oh — 



Cob 



t - 



cos 6* r 



k 



(6) 
(7) 



where k = g^^/l — (3q, and (3q is the initial velocity. 



/Xob can also be written in terms of tob directly. If we eliminate t in Equations 6 and 7, 
we obtain 



A*ob 



1 (y^^_T.), 



ksinO 



(8) 



where 



Ti — 1 — tohk cos 9 — Po cos 9 
T2 = {tl^k'' + 2Ukf3o)sm^9. 
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The description of the other parameters of the jet components, such as apparent velocity 
Papp and Doppler boosting factor D = 7~^(1 — cos^)"-*^, can also be derived from Equations 
6, 7 and from the basic equations in Section 2.1. The corresponding formulas are 

g ^ sme{kt + f3o) 

^kH"^ + 2kPot + 1 - cos e{kt + /3o) 



y/kH^ + 2kPot + 1 - cos 9{kt + /3o) ' 



D = _ (10) 



Compared to the non-accelerating relativistic jet models, the accelerating model has 
some new interesting features. When t is infinitesimal, the initial apparent velocity will 
be — /3osin^/(l — PqCOsO). As t increases, the apparent velocity will approach to 
its maximum value = cot(^/2) (Blandford & Konigl 1979). Concerning the Doppler 
factor D, it will first increase with t, and it will reach its maximum value D^^, where 
f3 = cos 6 and = D. If the acceleration goes on, then D will decrease, to values even lower 
than 1. Therefore, in an accelerating relativistic jet, both Doppler boosting and deboosting 
phenomena can be observed. 

The LASM model can be applied to estimate the parameters of a relativistic jet. Equa- 
tions 6 and 7 show that the shape of a proper motion curve is controlled by three parameters, 
i.e. (3q, g, and the viewing angle 6. By comparing the model to the observed VLBl data, 
these parameters can be calculated. If multi frequency flux density data were combined with 
proper motion data, then the model fitting results would be more accurate. 

Since the relativistic jets in AGNs probably exhibit different structures on different 
scales, we should carefully consider where the LASM model can be applied. On the large 
kpc scale, the jets usually have collinear structures. Although superluminal motions were 
detected by HST observations on such scale in M87 (Biretta et al. 1999), there are not enough 
data for a quantitive estimate of the jet parameters with LASM model. However, in this 
situation, the LASM model can still be used to qualitatively explain some interesting profiles 
of the jets (see Section 4.2). On the scale of tens of pc, the jets probably show 'wiggles' or 
hehcal structures, and the LASM model cannot be used there. On the small scale, especially 
in the core regions of jets, where jet components usually move out coUinearly, the LASM is 
applicable to assess the physical and geometrical parameters of jets. 



3. Application to 3C273 



3C 273 (J1229-I-0203) is a low optical polarization quasar (LPQ) with z — 0.158 (Strauss 
et al. 1992). It was the first object to display apparent superluminal motion on parsec scales 
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(Pearson et al. 1981). The radio structure of 3C 273 shows a well defined core-jet morphology 
from the mas scale up to the arc-second scale. The jet extends out to 20 arc-seconds, with 
the ridge line of emission showing a clear 'wiggle' (Davis et al. 1985). VLBI imaging at 
frequencies from 5 to 100 GHz show that the ridge line of the jet is curved on a scale 
from 0.05 to 25 mas, oscillating around the main orientation of the jet (Baath et al. 1991; 
Krichbaum et al. 1990). Zensus et al. (1990) suggest that the discrepancy in the position 
angles of some of the components observed at different frequencies may reflect a spectral 
index gradient across the jet. 

Since 1990, 3C273 has been monitored with VLBI from 15 to 86 GHz. For the com- 
ponents with enough data points at small (< 2mas) and large (> 2mas) core separations, 
quadratic flts to the radial motion r{t) describe the observations much better than linear 
fits (Krichbaum et al. 2001). This suggests that these components indeed have accelerating 
proper motions. Furthermore, in the inner region with radial distance less than 8 milli- 
arcsccond, these components seem to move outward coUinearly. Therefore, the LASM model 
is suitable to describe the observed data, and estimate the viewing angles and accelerations 
the components via least square fits. 

The components Cll, C12 and C13 in Krichbaum et al. (2001), which have enough good 

quality data, arc chosen for the model fitting. The location of a component at the first epoch 
is set as the origin of the coordinate system, with tob = and /iob = 0. The original proper 
motion data and their fitted lines arc shown in Figure 1. The corresponding parameters are 
listed in Table 1. Figure 2 shows the predicted evolution of the Doppler factors D of the 
components. 

Model fitting results show that the observed proper motion curves are sensitive to the 
viewing angle 6 and to the initial velocity f3Q. These two estimated parameters have relatively 
small errors. The acceleration factor k, however, has a larger error. The reason is that when 
a component accelerates to an ultra-relativistic velocity, its apparent velocity will remain 
almost constant, with /3app ~ cot(^^/2), thus lateral observed data will not dramatically 
contribute to the estimate of k. The values are lower than 1. This may due to the fact 
that the errors of the original proper motion data have been over-estimated. Table 1 and 
Figures 1 and 2, show that the Doppler factors of all the three components are decreasing 
in the last epoch. 
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4. Discussion 

4.1. Viewing angle 

The viewing angle plays an important role in the unification for radio-loud AGN. It is 
believed that all radio-loud AGNs have a similar physical structure: a supermassive black 
hole (probably spinning) located at their centre, whose gravitational potential energy is 
the ultimate origin of the AGN luminosity; an accretion disk surrounding the black hole, 
with relativistic jets emitted perpendicular to it; broad line and narrow line regions, and 
a dusty torus located at a larger distance from the innermost nucleus. The appearance of 
AGNs, however, depends so strongly on orientation that our current classification schemes 
are dominated by random pointing directions instead of more grounded physical properties. 
Thus, to understand the intrinsic mechanism of the AGNs, beaming effects induced by the 
viewing angle should be removed (Urry & Padovani 1995). 

Besides the apparent velocity measured by VLBl, other observational properties are 
needed to estimate the viewing angle. The apparent velocity /3app = (3 sin 6*/ (1 — /3 cos 9) de- 
pends both on 9 and (3, so 9 cannot be derived only from /9app- K both apparent velocity /3app 
and Doppler boosting factor D are known, then the viewing angle can be calculated (Ghis- 
ellini et al. 1993). /3app can be directly observed by VLBL Gurrently, there are three ways 
to estimate D. The first one is based on synchrotron self-Compton (SSC) model (Marscher 
1987). The second one depends on equipartition brightness temperature (Readhead 1994). 
The last one can be derived from the flux density variability of blazars (Lahteenmaki & 
Valtaoja 1999). Compared to the LASM model, all of these three models need extra physi- 
cal parameters which are probably difficult to assess. For example, SSC model requires the 
flux at the turnover frequency, which is usually substituted by the observed VLBI frequency 
in most cases, with an additional source of errors. Furthermore, these three models just 
provide the average values of the Doppler factors of jets, while the LASM model can trace 
the evolution of the Doppler factor as well as of the apparent velocity of the individual 
components. 

The LASM model provides a direct and rehable way to estimate the viewing angle of a 
relativistic jet, provided that there are enough high angular resolution VLBI monitoring data. 
As shown in Section 2.2, the shape of a proper motion curve is controlled by the viewing 
angle and other two parameters. Through nonlinear fitting, all the parameters could be 
evaluated. Furthermore, the acceleration will also influence the observed flux densities of jet 
components (see Equation 10). If flux density data are incorporated into model fitting, the 
accuracy of the estimation of viewing angle should be improved. 

The model fitting results of 3G 273 show that the jet components Gil G12 and G13 
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have different viewing angles (See Table 1). This is clear from the VLBI images (Krichbaum 
et al. 2001), where different jet components move out along projected trajectories with 
different position angles. The maximum difference of position angles can reach 40 degree. 
The difference of position angles as well as viewing angles may indicate that there probably 
exists a super-massive binary black hole system (Sundelius et al. 1997) where the orbital 
motion of the relativistic jet emerging from one massive black hole causes the variation of 
viewing angles (Rieger & Mannheim 2000). Another possible explanation is precession 
(Gower et al. 1982), and the precession angle and period can be used to estimate the masses 
of the black holes (Romero et al. 2000). 



4.2. Doppler factor 

During the accelerating stage, the evolution of the Doppler factor will affect the flux 
density variability of a jet component. The observed flux density Soh of a component is 
related to its intrinsic flux density S by S'ob = SD^'^"' (Blandford & Konigl 1979), where a 
is the spectral index. A small variation of D will cause a large change of the observed flux 
density. As seen in Figure 2, during the acceleration, the evolution of a Doppler factor is 
divided into two stages. In the first stage, D increases as the intrinsic speed of a component 
increases. In the second stage, as the intrinsic speed keeps raising, D will gradually decrease, 
going to values lower than 1. Assuming a = 0.5, the observed flux densities of Cll, C12 and 
C13 are boosted up to the maximum level of 1267, 140 and 820 times respectively. If we 
set D = 1 as the detection threshold of VLBI, then the traceable lengths of the components 
Cll, C12 and C13 range roughly from 3 mas to 60 mas. This is in agreement with the 
observations. 

These considerations allow us to interpret the gaps between the inner and outer jets of 
some radio loud AGNs, such as for instance 3C 273 and M87. For 3C 273, the inner jet is 
about 100 mas long at 1.67GHz (Davis et al. 1991). The outer jet is about 13" away from 
the inner jet and extends about 10" (Sambruna et al. 2001; Marshall et al. 2001). Inside 
this gap, no radiation is detected. The possible reason is that continued acceleration of 
components will cause Doppler dcboosting. Therefore, at some distance, these components 
will disappear under the detection level of present observation. However, the acceleration 
will end when the internal energy is mostly converted to kinetic energy. From then on, due 
to the interaction between the jets and the medium of their host galaxy at sub-kpc scales, 
the components will decelerate and and their Doppler factors will increase. Consequently, 
the Doppler boosted components will appear again in the outer jet region. 
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Fig. 1. — LASM model fitting results of the components Cll, C12 and C13 in 3C 273. 
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Fig. 2. — The evolution of the Doppler factors of the components Cll (Dotted hne), C12 
(Dashed hne) and C13 (FuU hne) in 3C 273. 
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Table 1. LASM model fitting results of 3C 273 with confidence level of 68.3%. DOF is the 
number of data used for fitting, 6 is the viewing angle of a component, is the initial 

velocity, and k is the acceleration factor. /3°pp and /3£pp, 7° and 7-'^, and denote the 
initial and final (at last observation epoch) apparent velocities, Lorcntz factors and 
Doppler boosting factors respectively. _D™"^ is the maximum Doppler factor where 

(3 = cos 6 and j = D. 



Parameters 


Cll 


C12 


C13 


DOF 


26 


26 


29 


6 (degree) 


7.5±0.3 


14.2±0.2 


8.4±0.2 




0.954±0.009 


0.7±0.2 


0.937±0.015 


k 


0.009±0.003 


0.2±0.1 


0.03±0.01 


/^app 


2.29 


0.53 


1.88 


11.2 


8.0 


12.2 


7° 


3.2 


1.4 


2.9 




12.7 


39.6 


20.3 




5.4 


2.2 


4.8 




6.8 


0.82 


4.1 


jymax 


7.7 


4.1 


6.8 




0.71 


0.34 


0.78 



